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Abstract 

Strong X-ray flares from the blazar Mrk 421 were detected in 2010 January and 
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February through the 7 month monitoring with the MAXI GSC. The maximum 2 
- 10 keV flux in the January and February flares was measured as 120 ± 10 mCrab 
and 164 ± 17 mCrab respectively; the latter is the highest among those reported from 
the object. A comparison of the MAXI and Swift BAT data suggests a convex X-ray 
spectrum with an approximated photon index of T > 2. This spectrum is consistent 
with a picture that MAXI is observing near the synchrotron peak frequency. The 
source exhibited a spectral variation during these flares, slightly different from those 
in the previous observations, in which the positive correlation between the flux and 
hardness was widely reported. By equating the halving decay timescale in the January 
flare, t& ~ 2.5 x 10 4 s, to the synchrotron cooling time, the magnetic field was evaluated 
as B ~4.5 x 10~ 2 G (5/10)" 1 / 3 . Assuming that the light crossing time of the emission 
region is shorter than the doubling rise time, t r ^2x 10 4 s, the region size was roughly 
estimated as R < 6 x 10 15 cm (5/10). These are consistent with the values previously 
reported. For the February flare, the rise time, t r < 1.3 x 10 5 s, gives a loose upper limit 
on the size as R < 4 x 10 16 cm (5/10), although the longer decay time ^ ~ 1.4 x 10 5 
s, indicates B ~ 1.5 X 10~ 2 G (5/10)" 1 / 3 , which is weaker than the previous results. 
This could be reconciled by invoking a scenario that this flare is a superposition of 
unresolved events with a shorter timescale. 

Key words: galaxies: BL Lacertae objects: individual (Mrk 421) — galaxies: 
active — X-rays: galaxies — radiation mechanisms: non-thermal 

1. Introduction 

Blazars, including BL Lacerate (BL Lac) objects, are a class of active galactic nuclei 
(AGN), from which a relativistic jet is emanating close to our line of sight. In addition to their 
non-thermal radiation ranging from the radio to gamma-ray frequencies and strong polarization 
at radio and optical frequencies, one of the most outstanding characteristics of blazars is their 
rapid and high-amplitude intensity variations or flares, which provide clues to the flow dynamics, 
as well as to particle acceleration and cooling processes operating in the jet. Since the Ginga 
(e.g., Sembay et al. 1993; Tashiro et al. 1995) and ASCA era (e.g., Takahashi et al. 2000; 
Tanihata et al. 2001), X-ray observations have been one of the most useful tools to study the 
variability of blazars. In addition, remarkable progress was successively accomplished by X-ray 
observations with RXTE and Swift observations (e.g., Kataoka et al. 2001; Tramacere et al. 
2009). 

Located at the redshift of z = 0.031, Mrk 421 is a high energy peaked BL Lac object 
(HBLs; Padovani & Giommi 1995). The object is one of the brightest extragalactic source at 
the very high energy (VHE) gamma-rays above ~ 100 GeV (Punch et al. 1992). A number of 
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studies (e.g., Kino et al. 2002) revealed that the mult i- frequency spectral energy distribution 
(SED) of the source and its variation are basically well understood within the framework of a 
simple one-zone synchrotron-self- Compton (SSC) model (Band & Grindlay 1985). For Mrk 421, 
the X-ray and VHE gamma-ray bands correspond to, or are located slightly above the peak 
frequencies of the synchrotron and inverse Compton (IC) spectral components, respectively. 
As a result, Mrk 421 is one of the most extensively studied blazars in the X-ray band (e.g., 
Takahashi et al. 2000; Tanihata et al. 2004), and X-ray flares were frequently reported (e.g., 
Tramacere et al. 2009; Donnarumma et al. 2009). 

As the first astronomical mission on the Exposed Facility of the Japanese Experiment 
Module "Kibo", attached to the International Space Station (ISS), Monitor of All-sky X-ray 
Image (MAXI; Matsuoka et al. 2009) started its operation in 2009 August. Thanks to its 
unprecedentedly high sensitivity as an all-sky X-ray monitor, and to its capability of real- 
time data transfer, MAXI is able not only to make a continuous monitor of X-ray sources 
including AGN, but also to alert quickly various transient X-ray phenomena, like flares of 
blazars. Actually, in the 7 month MAXI observation, we have successfully made quick alerts of 
two strong X-ray flares from Mrk 421 (Isobe et al. 2010a; Isobe et al. 2010b); in one of these 
flares the object exhibited the highest X-ray flux among those ever recorded. In the present 
paper, we report the X-ray features of Mrk 421 in these flares, and discuss their implication. 

2. Observation and Results 

Among the two X-ray instrument onboard MAXI, the Gas Slit Camera (GSC) and the 
Solid-state Slit Camera (Tsunemi et al. 2010), we analyze only the GSC data in the present 
paper, because of its higher sensitivity and sky coverage (~ 97% per day). The GSC utilizes 
12 position-sensitive large-area proportional counters sensitive to X-ray photons in 2 - 20 keV, 
although 4 out of them were switched off, on 2009 September 22, due to unexpected discharge 
events. 

The MAXI GSC signals from Mrk 421 were integrated orbit by orbit within a 3° x 3° 
square aligned to the scan direction centered on the source, while the background level was 
evaluated from two squares with offsets of ±3° toward the scan direction in the sky, of which 
the size is same as that of the source region. We normalized the background-subtracted signal 
counts by a time-integrated effective area of the collimator plus slits above the proportional 
counters. In the following, the data from all the activated counters are summed up. 

As assistance for the interpretation of the MAXI GSC data, we utilized the hard X- 
ray data taken with the Burst Alert Telescope (BAT; Barthelmy et al. 2005) onboard Swift. 
The Swift BAT is composed of the coded aperture mask and large area (5200 cm 2 ) CdZnTe 
detectors. Thanks to its wide field of view (~ 1.4 str), the typical sky coverage of the instrument 
exceeds ~ 70 % (Krimm 2009) in one day observation. The daily Swift BAT lightcurves of Mrk 
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421 in the 15 - 50 keV range are available from the Swift BAT transient monitor results x , 
provided by the Swift BAT team. 

Figure 1 (a) and (b) show the normalized signal count rates from Mrk 421 derived with 
the MAXI GSC in 2 - 4 keV and 4-10 keV, F 2 _ 4 and F 4 _i respectively, from 2009 August 22 
(MJD = 55065) to 2010 March 20 (MJD = 55275). On the plot, 10 % of the count rate from 
the Crab Nebula (100 mCrab, F2-4 = 0.13 counts cm" 2 s" 1 and F4-10 = 0.125 counts cm -2 s" 1 ; 
Nakahira et al. 2010) is indicated by the dotted lines. We notice some relatively long data gaps, 
such as on 2009 September 3-7 (MJD = 55077 - 55081), November 12 - 16 (MJD = 55147 - 
55151), and 2010 January 19 - 26 (MJD = 55215 - 55222). During these periods, Mrk 421 was 
located in the direction toward the ISS orbital pole unobservable with the MAXI GSC. The 
Swift BAT lightcurve in the 15 - 50 keV range -F15-50 is compared in figure 1 (c). 

The lightcurve reveals that Mrk 421 is highly variable in all the energy ranges. The 
fractional root mean square variability parameter (e.g., Rodriguez-Pascual et al. 1997) in F2-4, 
F 4 _ 10 and -F 15 _ 50 were evaluated as 0.52, 0.71 and 0.70, respectively, from the daily averaged 
lightcurves. The GSC has detected at least four flares with an X-ray flux of > 100 mCrab, on 
2009 November 1 (MJD = 55136), 2009 November 11 (MJD = 55146), 2010 January 1 (MJD = 
55197; Isobe et al. 2010a), and February 16 (MJD = 55243; Isobe et al. 2010b). Unfortunately, 
the observational conditions for the first and second flares were found to be relatively bad, 
although the second flare appears to be concurrent with the hard X-ray one, detected with 
the Swift BAT on 2009 November 10 (MJD = 55145; Krimm et al. 2009). Therefore, we 
focus on the GSC data associated the third and forth flares which are indicated by two arrows, 
denoted as Epoch A and Epoch B respectively in figure 1 (b). During Epoch A, no significant 
increase in -Fi5_5o was apparently found, while the object brightened in all the energy ranges 
during Epoch B. The flare of Epoch A (MJD = 55197) was not covered with the Swift BAT, 
due to sparse sampling toward Mrk 421 during this period. In addition, the hardness variation 
between Epoch A and Epoch B (before the first peak) was found to be statistically insignificant, 
as we show below. 

Figure 2 shows the 2-10 keV GSC lightcurve, F2-4 + -F4-10, during Epochs A and B, 
together with hardness ratios, defined as HR1 = -^4-10/^2-4 and HR2 = Fi 5 _ 5Q / F2-4. The 
lightcurve in Epoch A peaks at MJD = 55197.4, with a 2 - 10 keV flux of 120 ± 10 mCrab, 
averaged over 6 hours around the peak. Here and hereafter, the error represents the la statis- 
tical one. The doubling time scales in the rise and decay phases of this flare were estimated as 
t r ^2x 10 4 s and t& ~ 2.5 x 10 4 s, respectively. During Epoch B, the source activity continued 
for about one month, and we found several flares in the lightcurve. The peak 2-10 keV flux 
of the first flare in Epoch B was measured at MJD = 55243.6 as 164 ± 17 mCrab on 6 hour 
average. The decay time of this flare is derived as td ~ 1-4 x 10 5 s, a factor of 5 - 6 longer than 
that of the Epoch A flare. Due to the data gap with > 1 day just before the peak seen in figure 

1 http: / /heasarc. gsfc.nasa.gov / docs / swift / results /transients / 
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2, we put only a loose upper limit on the rise time as t T < 1.3 x 10 5 s. The lightcurve exhibited 
the second peak at MJD = 55249.6 with a 2 - 10 keV flux of 108 ± 7 mCrab. 

Utilising the daily averaged hardness ratio, we investigated quantitatively the spectral 
variation of Mrk 421 during Epochs A and B, as summarised in table 1. The hardness of the 
source was found to stay unchanged during the Epoch A within the statistical uncertainties 
(x 2 /dof = 7.3/18 and 14.9/15 for HR1 and HR2, respectively). The average and standard 
deviation of HR1 were derived as HR1 = 0.74 and <j\ = 0.10, respectively, while those of HR2 
were calculated as HR2 = 5.3 x 10 -2 , and cr 2 = 2.3 x 10~ 2 . Before the decay of the first flare in 
Epoch B (MJD = 55232 - 55245; corresponding to the arrow denoted as "Before" in figure 2), 
the object exhibited a hardness similar to that in Epoch A. No statistically significant spectral 
variation was observed during this phase (x 2 /dof = 6.6/11 and 7.4/9 for HR1 and HR2). 
After the decay of the Epoch B first flare (MJD = 55246 - 55275; indicated by the arrow 
with "After" in figure 2), the spectrum of the source seems to harden into HR1 = 0.90 and 
HR2 = 10.0 x 10~ 2 , with the standard deviation of <j\ = 0.15 and <Ji = 2.7 x 10~ 2 , respectively. 
In addition, the source hardness is found to be variable (x 2 /dof = 42.1/27 and 73.8/27 for HR1 
and HR2) during this period. Toward the end of Epoch B (MJD > 55264), the source appears 
to recover a hardness similar to those before the first peak (and hence those in Epoch A). In 
order to visualize these spectral variation associated with the flares, HR1 and HR2 are plotted 
against the source X-ray intensity, -F2-4 + -^4-10) m figures 3. 

Figure 4 shoes the color-color plot between HR1 and HR2 with a time resolution of 3 
days. From this figure, we roughly evaluated the shape of the X-ray spectrum of Mrk 421. If 
a simple power-law (PL) model was assumed, the observed values of HR1 correspond to a 2 
- 10 keV photon index of T = 2 - 2.5, as shown by the vertical dashed lines in figure 4. The 
systematic uncertainty in the photon index due to the current GSC response is estimated as 
Ar = 0.1 (Nakahira et al. 2010). The figure also indicates a spectral softening from the MAXI 
GSC to the Swift BAT energy ranges. 

3. Discussion 

As shown in figures 1 and 2, at least the two active phases of Mrk 421 in 2010 January 
(Epoch A) and February (Epoch B) were extensively monitored with the MAXI GSC. The 
maximum 2-10 keV flux of the flare in Epoch A was derived as 120 ± 10 mCrab at MJD 
= 55197.4. During Epoch B, the MAXI GSC revealed a significant long-term activity of the 
source for nearly a month with a 2 - 10 keV flux of ^ 50 mCrab exhibiting multiple flares. 
This demonstrates a great advantage of the MAXI GSC for monitoring long-term variations of 
blazars. Especially, the first flare of Epoch B, of which the peak flux was measured as 164 ± 17 
mCrab (MJD = 55243.6), is probably associated with a strong and variable activity in the 
VHE gamma-rays with the maximum VHE flux exceeding the 10 Crab level, detected on 2010 
February 17 (MJD = 55244) by the VERITAS Observatory (Rene et al. 2010). 
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Since it was identified as the first extragalactic VHE 7-ray source in 1990s (Punch et al. 
1992), Mrk 421 has been one of the most extensively studied blazars. As a result, a number 
of X-ray observations of the object were conducted with various X-ray observatories, including 
Ginga (Makino et al. 1992), ASCA (Takahashi et al. 1996; Takahashi et al. 2000; Tanihata 
et al. 2001), RXTE (e.g., Kataoka et al. 2001), XMM-Newton (e.g., Tramacere et al. 2007), 
Suzaku (Ushio et al. 2009) and so forth. Except for in significant X-ray flares, its 2-10 keV 
X-ray flux was reported to be typically below ~ 50 mCrab. Thanks to recent X-ray monitoring 
observations with Swift and the All Sky Monitor onboard RXTE, strong X-ray flares were 
successively reported from the object; these include the flare in 2006 June with the 2-10 
keV flux of ~ 85 mCrab (Tramacere et al. 2009), the flare in 2008 June with ~ 130 mCrab 
(Donnarumma et al. 2009) and the flare in 2009 December with ~ 100 mCrab (D'Ammando et 
al. 2009). It is important to note that the maximum X-ray flux in the first flare of Epoch B, 
measured with the MAXI GSC, is higher than any other previous results. Therefore, we have 
concluded that the MAXI GSC discovered the strongest X-ray from Mrk 421. 

The color-color plot in figure 4 suggests that the MAXI GSC spectrum of Mrk 421 is 
roughly approximated by a PL with a photon index of T = 2 - 2.5. A convex curvature in the 
MAXI GSC and Swift BAT energy range inferred from a comparison of HR1 to HR2 in figure 
4, together with the previous determination of the synchrotron peak energy in Mrk 421 around 
the observed photon energies of E p = 0.3 - 5 keV (e.g., Tanihata et al. 2004; Tramacere et al. 
2007; Ushio et al. 2009), implies a smoothly curved X-ray spectrum of the source. 

In the case of HBLs, where the synchrotron peak frequency is located at the X-ray 
band, the timescale of the X-ray flux variation is thought to be controlled mainly by interplay 
between the following timescales (Kataoka 2000); the acceleration timescale of high energy 
electrons r acc , the radiative cooling time r coo i, and the light crossing time of the emission region 
r crs- We, here, describe all the timescales in the observer's frame, after converting those in the 
source rest frame, r', as r = r'/5, in order to relate them directly to the observed timescales. 
For a significant flare to take place, a condition of r acc < r coo i is required. Assuming a flare is 
a single event, the rise timescale t r is expected to be determined by the longer of r acc and r crs , 
while the longer of r coo i and r crs is thought to dominate the decay timescale t&. The soft X-ray 
lag, discovered with ASCA from Mrk 421 in the decay phase of an X-ray flare (Takahashi et al. 
1996), was regarded to be produced by the energy dependence of t coo i. In contrast, symmetric 
X-ray lightcurves with t r ~ observed in day-scale flares from several HBLs, were interpreted 
to be shaped by r crs (Tanihata et al. 2001). 

Because the 2-10 keV MAXI GSC spectrum around the maximum of the flares was 
found to be consistent with T ~ 2, we regard that the MAXI GSC observed near the synchrotron 
peak frequencies. Therefore, we adopt the above argument on the flare timescale to evaluate 
the physical parameters of Mrk 421 in these flares. In the flare of Epoch A and the first flare of 
Epoch B, td was measured to be longer than t r . This indicates that t r and td are dominated by 
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the different timescales with each other. As a result, we speculate a condition of r coo i > r crs > r acc 
or T C ooi > r acc > r crs in these flares, corresponding to t d ~ r coo i and t r > r crs . 

The synchrotron radiation usually dominates other cooling process in HBLs, includ- 
ing Mrk 421 (e.g., Kataoka 2000). Assuming an isotropically distributed electrons in a ran- 
domly oriented magnetic field, the synchrotron cooling time scale is expressed as r coo i ~ 
1.5 x 10 3 fi~ 3 / 2 -E k ^y /2 ^~ 1//2 s (re-evaluated from Takahashi et al. 1996), with B, E ke y, and 5 being 
the magnetic field in G, the synchrotron photon energy in keV, and the beaming factor of the jet, 
respectively. We roughly evaluated a time-averaged magnetic field of B ~ 4.5 x 10~ 2 (5/10)~ 1 / 3 
G in the flare of Epoch A, from the relation of t& ~ t coo i. Here, we assume the typical syn- 
chrotron photon energy of i?keV ~ 4 keV, by averaging the X-ray spectrum of V ~ 2 in the 
MAXI GSC energy range (2 - 10 keV). This is barely consistent with the value determined 
from the one-zone SSC modeling to the observed SED of Mrk 421 in previous observations 
within the systematic uncertainty (B = 0.036 - 0.44 G; Kino et al. 2002). A spectral softening 
or soft lag, due to the energy dependence of r coo i, was unobservable within the MAXI GSC 
signal statistics. On the other hand, we can put an upper limit on the size R of the emission 
region, as R < cr crs <5/(l + z), where c is the speed of light. The condition of t T > r crs gives 
R < 6 x 10 15 (5/10) cm. This upper limit appears to be consistent with the values reported in 
the previous strong X-ray flares (e.g., Donnarumma et al. 2009). 

In the similar manner, we put a loose upper limit from t T on the region size of the Epoch 
B first flare as R < 4 x 10 16 (5/10) cm. The magnetic field, B ~ 1.5 x 10- 2 (J/10)- 1 / 3 G, is 
found to be weaker by a factor of ~ 3 than that in the Epoch A flare due to td longer by a 
factor of 5 - 6, while the spectral softening (from HR2 = (7.6 ± 0.9) x 10~ 2 at MJD = 55243 
to (5.1 ±0.7) x 10^ 2 at MJD = 55245), suggested in figures 2 and 3 during the decay phase 
of the Epoch B first flare, appears to be related to the cooling process. We cannot deny a 
simple explanation for the possible weak magnetic field that this flare was produced in a part 
of the jet different from that of the Epoch A flare with slightly different physical parameters 
and/or conditions, from the MAXI data alone. However, we regard that td is not actually 
determined by r coo i in the case of this flare, because the derived magnetic field is found to be 
out of range of the typical value of Mrk 421 in the previous studies (e.g., Kino et al. 2002). 
Based on close examination on X-ray lightcurve of HBLs obtained in the long-look ASCA and 
RXTE observation with a duration of ^ 1 week, Tanihata et al. (2001) proposed a possibility 
that day-scale flares of HBLs are composed of short events with a timescale of ~ 10 4 s. A 
similar scenario is possible to account for the duration of the first flare of Epoch B (t& ~ 1.6 
day), longer by a factor of 5 - 6 than that of the Epoch A flare. It is important to note that 
in order to validate the physical quantities estimated here, examination of a multi-wavelength 
SED of Mrk 421, widely covering the synchrotron and IC components, and its variation are 
required, although it is beyond the scope of the present paper. 

Within the framework of a simple one-zone SSC model, the peak energy and luminosity 
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of the synchrotron component is thought to scale as E p oc B^ 2 5 and L p oc B 2 ^ 2-1 *^ 5 2+Tp N e V , 
where 7 P , T p (< 2), N e and V is the maximum Lorentz factor of the radiating electrons, the 
synchrotron photon index below E p (assuming a simple PL spectrum), the normalization of the 
electron number density spectrum, and the volume of the emission region, respectively. From 
these scaling relation, a change in B, 7 P or 5 is expected to be observed as a positive correlation 
between E p and L p , while only L p is anticipated to increase when N c or V (corresponding to the 
total electron number) increase. The clear positive E p -L p correlation, widely reported from 
Mrk 421 in the previous observations (e.g., Tanihata et al. 2004; Tramacere et al. 2007), is 
thought to suggest the change in B, 7 P or 5. 

The rise of E p is thought to be recognised as a rise in HR1 and HR2. The hardness- 
intensity relation, shown in figure 3, suggests a slightly complicated spectral variation, which 
appears not necessarily to follow the E p -L p positive correlation, especially in Epoch B. We 
notice that HR2 become larger, after the end of the first flare in Epoch B, as clearly shown 
in figure 2 and table 1. Although at the second flare of Epoch B the source is fainter by 
~ 30 % than at the first one, its hardness (HR2 = (10.3 ± 1.5) x 10~ 2 at MJD = 55249) is 
higher than that at the first flare (HR2 = (7.6 ± 1.5) x 10~ 2 at MJD = 55243). In addition, 
the hardest X-ray spectrum was realized, when the source flux was about one third of the 
maximum at the first peak (HR2 = (17.4 ±2.5) x 10~ 2 at MJD = 55262). In order to interpret 
such a spectral behaviour within a simple one-zone SSC framework, fine tuning of the physical 
parameters is required (e.g., a combination of the increase in B, 7 p or 5 and the decrease in 
N c or V). Otherwise, a multiple-zone/component SSC model (e.g., Ushio et al. 2009) could be 
preferred. For detailed examination, it is necessary to investigate the mult i- frequency SED and 
its variation. 
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Table 1. Summary of spectral variation of Mrk 421 



* Before the first flare (MJD 
t After the first flare (MJD = 

* Average of HRl and HR2. 
§ Standard deviation of HRl 



Epoch 


A 


B * 


B f 


HRl * 


0.74 


0.73 


0.90 


<T i § 


0.10 


0.13 


0.15 


X 2 /dof 


7.3/18 


6.6/11 


42.1/27 


HR2 i 


5.3 x 10~ 2 


6.2 x 10~ 2 


10.0 x 10~ 2 


er 2 § 


2.3 x 10~ 2 
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Fig. 1. Daily lightcurve of Mrk 421. The panels (a) and (b) show the normalized count rate in the 2-4 
keV (i<4_io) and 4-10 keV (-F2-4) ranges, respectively, measured with the MAXI GSC from 2009 August 
22 (MJD=55062) to 2010 March 20 (MJD=55275). For MJD = 55197, 55243 and 55249, F4-10 and F 2 -a 
are plotted in a time resolution of 6 hours, in order to resolve the peak. The Swift BAT lightcurve in 15 
- 50 keV, Fi5_5o, is plotted in the panel (c). All the errors represent the ler statistical one. The dotted 
lines indicate the count rate of 100 mCrab in the individual energy ranges. The two arrows in the center 
panel define Epochs A and B. 
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Fig. 2. Lightcurve of Mrk 421 during Epochs A (panel a) and B (panel b). The GSC 2-10 keV 
count rate F4-10 + -F2-4 in a time resolution of 6 hours is shown in the top. A time resolution of 1.5 
hour is adopted, when the 6 hour averaged flux is higher than the 100 mCrab level (the dotted line). 
The center and bottom display the daily-averaged hardness ratios, defined as HR1 = F4—10/F2-4 and 
HR2 = ^15-50/^2-4, respectively. For MJD = 55197 in panel (a), HR1 is displayed every 6 hours. The 
hardness for the Crab-like spectrum with the photon index of T = 2.08 (Toor & Seward 1974) with the 
Galactic absorption column density of JVh = 3.2 x 21 cm -2 (Kalberla et al. 2005) is expressed by the dotted 
lines. The long and short arrows in the center of panel (b) define the phases before and after the decay of 
the Epoch B first peak (MJD = 55246), adopted in figures 3 and 4. 
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Fig. 3. HR1 (a) and HR2 (b) plotted against ^2-4 + ^4-101 in a time resolution of a day. Only the 
data with F2-4 + i<4_io > 0.05 counts cm" 2 s _1 arc displayed. The error bars indicate the la statistical 
errors. The red circles show the data during Epoch A, while the green triangle and blue crosses represent 
the Epoch B data before and after MJD = 55246, respectively (indicated the two arrows in the center of 
figure 2). The data points covering the peaks of the Epoch A flare, the first and second flares of Epoch B 
are respectively indicated with Al, Bl and B2, respectively. 
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Fig. 4. Color-color plot between HR1 and HR2, averaged over 3 days. The symbol notation same as 
figure 3 is adopted for the individual phases. Those covering the peaks of the Epoch A flare, the first and 
second flares of Epoch B are indicated with Al, Bl and B2, respectively. The error bars indicate the la 
statistical errors. The dashed lines show the values of HR1 and HR2, for the photon index of T = 2.0, 2.5 
and 3.0, assuming that the Galactic absorption column density of TVh = 1-9 x 20 cm~ 2 toward Mrk 421 
(Kalberla et al. 2005). If the spectrum is a straight PL in the 2-50 keV range, the data points lie on the 
solid line. 
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